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Abstract 

We take graphene as an example to demonstrate that the present widely adopted expression is 
only the scattering component of a true 2D group delay in the condensed matter context, in which 
the spatial Goos-Hanchen (GH) shift along an interface contributes an intrinsic component. We 
relate the dwell time to spin precession and derive a relation between the 2D group delay and dwell 
time, whereby we for the first time reveal that, the group delay for 2D ballistic electronic systems 
can be directly observed by measuring a conductance difference in a weak-field spin precession 
experiment. This physical observable not only implies the group delay being a relevant quantity 
even in the condensed matter context, but also provides an experimental evidence for the intrinsic 
effect of the GH shift. Finally, we revisit the 2D Hartman effect, a central issue of the group delay, 
by analytically solving it via the vested relation and calculating the proposed observable at the 
Dirac point. 
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I. INTRODUCTION 



The tunneling of a particle through a barrier is one of the most ubiquitous and fundamen- 
tal quantum processes. Over eighty years ago, it was suggested that there is a time duration 
associated with such a process 1|. Among the various proposed expressions 2|, y], Wigner- 



Smith delay or group delay {r g 



6] are two well-established times, 



,|5| and dwell time (r^) [5, 
which describe the reflection or transition of a pulse peak and the momentary capture and 
release of a tunneling particle, respectively. ID tunneling group delay stems from the phase 
shift introduced by scattering at the interfaces [t| • It can be given by the eigenvalues of the 
Wigner-Smith time-delay matrix, or equivalently expressed by the change rate of the phase 
shift with respect to angular frequency 



T. 



St 



~dE' 



(1) 



where r = |r|e^ r (t = |t|e J<?it ) and r 



■k x l, see, Fig. lb) are the reflection (transmis- 



sion) coefficient and corresponding phase shift, respective 



y. In the 2D tunneling case, there 
of the reflected or transmitted 



generally is a nonzero lateral Goos-Hanchen (GH) shift 

beam along corresponding interface ^l^] (see Fig. lb). Obviously, it intrinsically modifies 
the phase shift at the interface, thus plays an important role in the 2D group delay. This 
fact has been noticed in optics L4J but long-term ignored in the condensed matter context. 
Instead, Eq. (1) contai ning only the scattering component has been widely used in earlier 
reports (see, e.g., Refs. 15I-I22]). Accordingly, the 2D Hartman effect 23[, a central issue of 
the group delay, needs to be revisited. 

On the other hand, measurements of the tunneling process are essential for all attosecond 
experiments. Very recently, group delay measurements are eventually achieved in atoms and 



molecules by approaches with rich skills 



24 



25| . However, these approaches are not feasible 



in condensed matters due to the fact that the numerous electrons on the Fermi surface are 
hard to be detected individually. For a long term, this dilemma has made the group delay 
being considered to be of little physical significance in this context So, to propose a 
feasible approach in such a context is a fundamental subject and still a key challenge in 



the field. Very recently, the authors in Ref 
a Larmor clock 



26 



16] made an attempt in graphene by using 



27] . However, the analysis was restricted to the Dirac point and the 



improper definition of 2D group delay was adopted. 



In this work, we demonstrate the effect of the GH shift on the 2D group delay and propose 
an approach for measuring group delay in 2D ballistic electronic systems. We take graphene 
as an example and show that the GH shift contributes an intrinsic component which 



adds an asymmetric feature to the 2D group delay. We find that the 2D dwell time can 
be related to a weak-field spin precession experiment and we derive a relation between the 
group delay and dwell time, whereby we reveal that, the 2D group delay (and hence the 
intrinsic effect of the GH shift on it) can be directly observed through simple conductance 
measurements in this experiment. Such an approach provides a general tool for group delay 
measurement in 2D ballistic electronic systems, and is much easier to be realized than the 



ones used in Refs. 



24( and 25|. We at last analytically resolve the 2D Hartman effect in 



graphene and investigate the proposed observable at the Dirac point. 



II. INTRINSIC EFFECT OF GH SHIFT ON 2D GROUP DELAY 



Let us consider a 2D quantum tunneling through a barrier in graphene [29(] (see, Fig. 1). 
The electron possesses Fermi energy E and a central incident angle of a. In the stationary 
state description, the electron can be represented as a wave packet as a weighted superpo- 
sition of plane waves (each being a solution of Dirac's equation). The appearing locus and 
time of the packet peak (equivalently, the electron) are determined by the condition that 
the gradient of the total phase in the wave vector (k = k x x + k y y) space must vanish. This 



is similar to the optical case 30j. A comparison of the conditions in direction for the re 
fleeted (transmitted) and the incident beam provides corresponding lateral shift jlO, 12, 13], 



af = —d(j) r /dky =F l/2/c x and af = —d(p t /dk y [see Fig. 1(b)]. Here £ ±' stands for the A and 
B component of the graphene spinor, respectively. A comparison of the conditions in magni- 
tude gives the average 2D group delay in reflection and transmission, r| = hd((p^ + k y a^)/dE, 
where average values = —d<j)^/dk y (ay = o t for a symmetric barrier) have be used since 
the two components are simply added in the group delay. For asymmetric barriers, there is 
a difference between r r g and r*, and the 2D group delay is defined as r g = |£| 2r |- 111 view 
of Eq. (1), it can be rewritten as 

r 9 = rf + if*, (2) 

where is due to the scattering and r g H = |£| 2(J £ sma/v F results from the GH shift. It 
is noted that the GH component has been widely ignored in the condensed matter context, 



see, e.g., Refs. 15|-|22|. At time hd(j) t /dE the wave front of the transmitted beam reaches 



(I, 0) and then propagates freely with a velocity of Vp to the final position (/, a t ), which will 
cost a duration of at sin a jvp since the wave front is perpendicular to the propagation direc- 
tion. This picture also holds for the reflected beam. A weighted average of the transmission 
and reflection gives t^ h . 

Fig. 2 shows clearly the contribution of the GH shift to the 2D group delay. As is 
seen, the scattering component is symmetric about the center of the transmission gap (TG) 
(E/V = cos -2 a) due to the symmetry of 0^ about the center. Contrarily, the GH component 
is asymmetric about the TG's center, which stems from the behavior of the quantum GH 
shift. The GH shift has the same trend as the classical shift (as = /tan/3 with sin /3 = 
hvpk y /(E — V) the refracted angle) predicted by the Snell's law 12J, which is negative 
(positive) in the low (high) energy region (see Fig. 2). Accordingly, through the intrinsic 
term in the phase shift (<j)^ + k y a^), the GH shift not only quantitatively contributes a part of 
order of as sin a/ 'vp in the 2D group delay but also, qualitatively, results in the remarkable 
asymmetric feature of the 2D group delay. 



III. MEASUREMENT OF THE GROUP DELAY BY SPIN PRECESSION 



We now seek physical observable for the 2D group delay in graphene. The Larmor 



precession of t 
dynamics 16 



re electron spin in a magnetic field provides a clock for studying the electron 



26 



3l| . Here, we calculate the transmission probability (T = |t| 2 ) and in 



a geometry where the magnetic field is applied in the graphene plane [16J. (Note the 2D 
dwell time has the same definition as the ID one.) An explicit equality between the dwell 
time and an average transmission probability < T >= (T zy — T zy )/(T zy + T zy ) is found in the 
weak-field limit: UsTd =< T > \b->o (see, Fig. 3). Here ojb = g^sB is the Larmor frequency, 
g the gyromagnetic factor, /ig the Bohr magneton, and B the magnetic field. The weakness 
of the magnetic field can be described by its reduced strength, B = Jtwb/2E . T z ^) y is 
the transmission probability for an electron incident with ^/-directed spin and transmitted 
with z(2)-directed spin. This equation holds for all Fermi energies and incident angles 
and can be interpreted physically in following. Multiplying the right side by H/2 gives the 
expectation value of S z , which is just determined by the Larmor frequency and the time 
the precession persists (obviously, the dwell time rather than r^* shown in Ref. [if]]) as 



the left side expressed. Note T implies an experimental observable, i.e., a conductance 
G(E) = G f^ 2 T(E,a)d(sma), with the conductance quantum of G = 2e 2 /h accounting 
for a twofold valley degeneracy. Thus, the above equality can be rewritten as 

T d (E,a)T(E,a)d(sma) = ~ ^v( E ) ^ (3) 

J-n/2 ^B^O 

which clearly indicates that the dwell times can be related to a spin precession experiment. 

We now try to relate the group delay to the conductance by obtaining the relation between 
it and the dwell time. Due to the 2D feature of the tunneling and spinor nature of graphene, 
the derivation for the relation is much more skillful and complex than the one for ID case of 
normal Fermions shown in Ref. Q. The detailed derivation can be found in the Appendix, 
here we give the concise result for a rectangular potential barrier 

Td = T g + Ti, (4) 

where r, = h\Re(r) cos a + Im(r) sin a] sina/E'cos 2 a is a self-interference delay stemming 
from the interference of the incident and reflection wave functions [2l, l32|. Note this term has 



so far been widely believed to disappear 17H20|. where the GH component was ignored and 
an improper variation method that is simply extended from Ref. j3] was taken. Fig. 4 shows 
the group delay, dwell time, and self-interference delay in reduced form as a function of the 
Fermi energy at a fixed incident angle, where e = E/E is a reduction factor. The reduced 
self- interference delay (— er») achieves the maximum at the TG's center and oscillates around 
zero outside the gap. The self-interference delay itself is important only in the low energy 
region (diverging as E" 1 when E — >■ 0), and disappears at (anti)resonant tunneling since 
there is no interference in front of the barrier. Accordingly, the group delay nearly coincides 
with the dwell time except within the low energy ranges or around the TG's center. 

The average self-interference delay, TjT<i(sina) oscillates with E and disappears 

at a relatively high Fermi energy (see the insert in Fig. 4). Therefore, Eqs. (3) and (4) 
imply that, for any Fermi energy compared with V, the 2D group delay can now be directly 
observed by spin associated conductance difference in a spin precession experiment 

™7r/2 /~< /~i 

-it/2 

As is seen in Fig. 5, the spin associated conductance difference which can be obtained by 



^ 2 r g (E,a)T(E,a)d(sina) « % -fe | B _> . (5) 



measuring G zy and G 2y (see insert in Fig. 5) increases for weaker B. At B = 10 3 , it is 



already a rather good measurement of the 2D group delay for E > 0.2V. Thus far, we have 
provided a feasible and relatively easy approach for measuring group delay in graphene, 
which can be applied to other 2D ballistic electronic systems, such as normal semiconduc- 
tors, bilayer and multilayer graphene, surface states of 3D topological insulators, and their 
heterojunctions. (Superconductors will be an exception as superconducting states may be 
destroyed by the magnetic fields). Since the experimental measurement is feasible, the group 
delay in ballistic electronic systems should be regarded as a relevant quantity. Meanwhile, 
the scattering component is larger (smaller) than the expected value of the group delay when 
E < V (E > V), a result of the GH shift of different signs. Then a comparison between 
the observed value of the conductance difference (right side of Eq. (5)) and the theoretical 
prediction of the average group delay (left side of Eq. (5)) can be utilized to probe the 
intrinsic effect of the GH shift. 



IV. THE 2D HARTMAN EFFECT IN GRAPHENE: ANALYTICAL RESULTS 
AND PHYSICAL OBSERVABLE 

We now revisit the 2D Hartman effect in graphene analytically. We find that the self- 
interference delay can be related to the difference in the expectation values of p x at the 
two outer or inner boundaries of the barrier, Tj = —h 2 {ip*dip j 'dx)\ l 011 /2j in p 2 . This can be 
rewritten as n = q 2 x T 2 (i)/p 2 x for q 2 x > (q 2 x < 0), where ri (2 ) = Pi( 2 )/jin, jin = v F cosa is the 
flux of incident particles, and -Pi(2) is the intensity (interference) component of the stored 
probability. Based on this equation and Eq. (4) (note = T\ + r 2 ) the group delay can be 
rewritten as r g = 7~i( 2 ) + (1 — A 2 )r 2 (i)/ cos 2 a for q 2 > (q 2 < 0), where A = (E — V)/E is a 
ratio of the kinetic energies inside and outside the barrier (see the appendix). The probability 
density for the evanescent mode (hi 2 = —q 2 > 0) simply decays exponentially. Taking 
the analytical expression of -Pi(2) by the transfer matrix method, one can find in the limit 
I —> oo, Pi — y E cos 2 a/Vn, and P 2 — * 0. Thus the dwell time, the group delay, and the self- 
interference delay are all saturated in this exponential limit: — > (nvp)~ l E cos a/V, r g — > 
(kwf)~ 1 (1 + A) / cos a, and Tj — > (nvp)' 1 ^ cos 2 a — (1 + A)F]/V cos a. The analytical results 
obtained here clearly indicate that the 2D Hartman effect does happen in graphene provided 
q x becomes imaginary under the barrier, which is the same condition for semiconductor-based 
structures with parabolic dispersions [7| • This result is not only quantitatively different from 
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Ref. (17| where only the scattering component was considered, but also essentially opposite 
to the Dragomans' conclusion that the Hartman effect is absent in graphene as declaimed 
in Ref. Q- 

We at last devote ourselves to the barrier-length dependent behavior of the proposed 
physical observable. Since there are no transmission modes, the average group delay is also 
saturated with the barrier length at the Dirac point (E = V) (see, insert of Fig. 6). This 
behavior and the saturation value can be directly detected by conductance measurements. 
Moreover, the group delay is dominated by the scattering component, which can be well 
approximated by Eq tanhL/i?sina with L = El sin a/ Eq. As can be seen in Fig. 6, al- 
though the GH component dramatically increases as \a\ and finally dominates the group 
delay, it makes negligible contribution to the average group delay, since the transmission 
probability (which can be well approximated by 1/ cosh 2 L) exponentially decreases with 
\a\. In other word, the intrinsic effect of the GH shift should be probed deviating from the 
Dirac point, which coincides with Fig. 5. In view of Eq. (5), the saturation value is given 
by (Eo/E) J^tanh L/(L cosh 2 L)dL = U((3)E /(n 2 E). For our case, E/Eq — 37r gives a 
value of about 0.18, which coincides well with the numerical result in Fig. 6. 



V. CONCLUSIONS AND REMARKS 



In summary, we have demonstrated the intrinsic contribution of the GH shift to the 2D 
group delay in graphene. More importantly, we have proposed a feasible and relatively 
easy approach for the observation of group delay (and hence the intrinsic effect of the 
GH shift) by conductance measurements in a weak-field spin precession experiment. We 
have also analytically revisited the 2D Hartman effect, and investigated the barrier- length 
dependent behavior of the proposed observable. Although quantitative results may vary for 
other systems which differ in electronic elementary excitations, the intrinsic effect of GH 
shift generally hold for any 2D ballistic electronic systems, and the proposed conductance 
difference provides a universal physical observable for 2D group delay in these systems 
(except superconductors). Besides, our results also imply that to construct a generalized 
Wigner-Smith time-delay matrix that is valid for 2D case and 3D case is an urgent issue in 
the condensed matter physics field. 
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Appendix A: Derivation for the relation between 2D r g and in graphene 

Let us begin with the single-particle Dirac equation that governs the low-energy excitation 
in graphene. In the barrier region it reads as [vpcr ■ p + V(x)]Sfr(x, y) = E^(x,y), where 
the pseudospin matrix cr has components given by Pauli's matrices and p = (p x ,Py) is 
the momentum operator. The eigenstates &(x, y) are two-component spinors with each 
component being the envelope function at sublattice site A/B of the graphene sheet. Due 
to the translational invariance along the y-axis, the wave function can be separated as 
VP = \iPa(x), ipB{x)} T e lkvy with k y = E sma/hvp. The two x-part components are related 
by a pair of coupling first-order equations 

d E - V 

— i) A = k y ^p A — — ij} B , (Ala) 

d E-V 

*q IpB = -ky^B TT VA, ( Al b) 

ox invp 
which implies a decoupled second-order equation for both the A-site and -B-site components 

d 2 



+ (j^) 2 l(E - Vf - E 2 sin 2 a] J ^ = 0, 



(A2) 



v dx 2 
where ^ = i/j a ,b- 

We carry out the energy-variational form and conjugate form of Eq. (A2) and upon 
integration over the length of the barrier we find 

(dto d^* _ O 2 ^ \ x=l _ f 2E C os 2 a-2V{x) 

{dE~dx- dEdx-J x=0 ~ J (hv F f ^ ^ (A3) 

It is seen that the left (right) part can be related to the group delay (dwell time), when is 
evaluated by the wave function outside (inside) the barrier. Note Eq. (A3) is only valid 
inside the barrier, we express the position derivation of the component inside the barrier by 
Eq. (Al) and their conjugate form, which then can be replaced by the corresponding ones 
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outside the barrier. Then the left part of Eq. (A3) can be evaluated. For the A-component 
it reads as G + F + cosct(— ir + ir*)/hvp, and for the B- component, the result becomes 
G + F + cosa(ie- i2a r-ie i2a r*)/Hv F , where G = ^{[B(0)-A(0)]\r\ 2 (j)' r + [B(l)-A(l)}\t\ 2 (t)' t }, 
F = JjL{[B(0) - A(0)]\r\\r\' +[B(l) - A(l)]\t\\t\'}, and the relation of lossless barriers \t\ 2 + 
|r| 2 = 1 is used. The notations are adopted as O' = dO/dE, A(x) = sin a + iX(x)e ia , 
B{x) = sin a — iX(x)e~ ta , and X(x) = E ~ V E ^ > , a ratio of the kinetic energy inside and outside 
the barrier. Since ip*ifj = if)\ipA + ^b^b, the relation for each component should be added 
to get the relation for the two-component spinor, which at last reads 

/q[A(x) — sin 2 a}\ifj(x)\ 2 dx 
Vp cos a 

=A(0)|r| 2 ^ + A(/)|t| 2 ^ - i\(0)H\r\\r\'- ( A4 ) 

. w , s ., i, ^ fRe(r) cos a + Im(r) sinal sin a dk 
i\{l)h\t\\t\' + h [ -^ ^ J — . 

This is a general result relating the integral weighted probability density inside the barrier 
(left part) and the weighted energy- variational behavior outside the barrier (right part). As 
is seen, the factor A(x) has a critical role in the relation. 

To clearly relate the general result Eq. (lA4p with both and , we consider a restricted 
condition that X(x) is a constant under the barrier (i.e., a rectangular barrier). Note this 
condition is not necessary for normal semiconductors case a reflection of the spinor 
nature of graphene. In this case, the third and fourth terms on the right side of Eq. (1A4j) 
disappear due to the lossless condition of the barrier \t\\t\' + \r\\r\' = 0, and Eq. flA4j) can 
be rewritten in terms of and sub-relation reads 

Td(X — sin 2 a) = r^*A + Tj cos 2 a, (A5) 

where a self-interference delay is found from the last term of Eq. (A4), n = h[Ke(r) cos a + 
Im(r) sin a] sina/E'cos 2 a. It is seen that even the GH component is ignored, the "group 



17|-|20| 



delay" in graphene does not equal the dwell time as indicated in Refs. 

The angle-variation of Eq. (A2) is straightforward following a similar way and the vari- 
ation result gives the sub-relation between and t^ h under the same restricted condition 
of constant A. The result reads as 

Trfsin 2 a = TgX + Ti(X — cos 2 a), (A6) 
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which, to our best acknowledge, has not been noticed up to now. Making a simple addition 
of the sub- relations in Eqs. (A5) and (A6) and taking into account Eq. (2), we finally get 



Td = Tn 



(A7) 



For the normal incident or ID tunneling case (a = 0) the self-interference delay vanishes 
[see the factor sin a], since no reflected portion thus no interference happens in front of the 



barrier due to the Klein paradox 



33|]. 



The relation revealed in Eq. (A7) and the expression for the self-interference de 



ay also 



hold for the tunneling of massless Dirac particles in the topological surface states [34| where 
the real electron spin rather than the sublattice structure in graphene provides the Dirac 
structure. It is also an analogue of electromagnetic waves in the photonic band gap structure, 
of which the normal incident case was investigated in Refs. 
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(a) ' (b) 

FIG. 1. Sectional (a) and top (b) view of schematic diagrams for a particle quantum tunnel- 
ing through a potential barrier in graphene. The red solid and blue dashed lines stand for the 
trajectories of the A and B components, respectively. 
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FIG. 2. The 2D group delay and its scattering and GH components (in units of the equal time 
To = 1/vf) as a function of the Fermi energy at a = 20°. The magenta thin solid curve stands for 
the time associated with the Snell shift rg = erg sina/vF- The parameters of the potential barrier 
(which also hold in Figs. 2-5) are 1/1$ = 1 and V/Eq = 3tt with Iq and Eq = Hvf/Iq being a length 
unit and energy unit, respectively. 
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FIG. 3. The dwell time and S-dependent < T > /u^ asa function of the Fermi energy at a = 20°. 
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FIG. 4. The reduced group delay, dwell time, self-interference delay, and transmission probability 
as a function of the Fermi energy at a = 20°. Insert: the average self-interference delay versus the 
Fermi energy. 
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FIG. 5. The average group delay, its scattering component, and the independent conductance 
difference (i.e., the right hand of Eq. (5)) as a function of the Fermi energy. Insert: G zy (solid) 
and Gzy (dashed) could be directly measured in the experiment for B = 10 -1 . 
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FIG. 6. The group delay, its two components, and transmission probability as a function of a 
at the Dirac point of a barrier of V/Eq = 3ir and I/Iq = 1.5. Insert: the conductance difference, 
average group delay, and its scattering component versus the barrier length at the Dirac point for 
V/E = 3vr. 
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